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Abstract 

In this paper we treat the 2-level system interacting with external fields without 
the rotating wave approximation and construct some approximate solutions in the 
strong coupling regime. 

In this paper we consider the 2-level system interacting with external fields (a radiation 

field for example ) and treat this system without the rotating wave approximation (RWA). 

Usually we assume the RWA and solve the equation. This approximation is appropriate 

in the weak coupling regime, however in the strong coupling one it is not applicable. See 

for example [1 and [2j. How can we treat the system in the strong coupling regime ? 

Shahriar et al treated this problem in |3j, [I] and [H] and constructed some approximate 

solutions which contained what they called the Bloch-Siegert oscillation. See also , [Z] , 

jHj, as another approachs. 
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By the way we have treated some important models coming from quantum optics in the 
strong coupling regime, see JUj, P~2] and ^3] (also ^2] for an experiment in the strong 
coupling regime). See [HJ as a general introduction in this field. This method is also 
applicable in the case, so we construct some approximate solutions which are completely 
different from those of Shahriar et al. We believe that our solutions are more general. 

Let {a\, 02, 03} be Pauli matrices : 

1 

0-1=1 |,0- 2 

1 

and a + = (1/2) (cri + iai), 0- = (1/2) (01 — ia\). 

We consider an atom with two energy levels which interacts with external field with 
gcos(ut). The Hamiltonian in the dipole approximation is given by 

H = -y0 3 + g cos(wt)<n, cos(cjt) = (l/2)(e^* + e"^'), (2) 

where u is the frequency of the external field, A the energy difference between two levels 
of the atom, g the coupling constant between the external field and the atom. In the 
following we cannot assume the rotating wave approximation (which neglects the fast 
oscillating terms), namely the Hamiltonian given by 

H = -y^s + § (e ia V+ + e—V_) . (3) 

Mysteriously enough we cannot solve this simple model completely (maybe non-integrable) , 
nevertheless we have found this model has a very rich structure. 

We would like to solve the Schrodinger equation with the Hamiltonian (J2J) in the 
strong coupling regime g ^> A. Therefore we change in (J2J a role of the kinetic term 
and the interaction term like 

H = Hq- = g COs(u;t)0i - y0 3 . 

First we solve the interaction term which is an easy task. Before that let us make some 
preliminaries. 



Let Wbea Walsh-Hadamard matrix 

1/11 

W 



(4) 



V^l -1 

then we can diagonalize <j\ by using this W as 0\ = WdsW' 1 . The eigenvalues of <j\ is 
{1,-1} with eigenvectors 



|1> 



V2 
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V2 



We note that cr 3 1 1) = |— 1), a 3 \ — 1) = |1) and 

/ 



111 

|1><1| = « 

1 1 1 1 



1/1-1 

-i><-i| = - 

^1-11 



V 
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W = W(\0)(0\)W, 



W = W(\1)(1\)W, 



where |0) = (1,0)* and |1) = {0,1^. 
Then it is easy to see 

H = gcos(ujt)\l)(l\-gcos(ujt)\-l}(-l\, 
so the Schrodinger equation i-^ip = H ip can be solved to be 

^(t) = {e"^ sin ^|l)(l| +e l ^ sm{ut) \-l)(-l\}^o, 



where ipo is a constant vector. Since we want to solve the full Schrodinger equation 

d , ( . A 



#0 - ) V>> 



(5) 



(6) 



(7) 



we appeal to the method of constant variation. Namely under ip Q — > ipo(t) in © we 
substitute (|BJ) into ((Zj), then we obtain 



d , A 

l — tyn = X 
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— {e 2i - sm(aJ ' } 1 1) (- 1 1 + e" 2i - s ™ (art) | - 1) ( 1 1 } 
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so we have only to solve the equation 

A 



d 



--{ 

9. I 



Jli— sin(u>t) _ i „—2i—sin(u>t) r ^ 



where |0)(1| = a + and |1)(0| = a_. 

Here we note that the solution ip in (JZ) is given as follows 

^(t) = ^{e _< S" n(wt) |0)(0| +e i S siri(a; * ) |l)(l|}0(t) 

with in ((H). 

To solve the equation (JHJ) we set : 



(9) 



(10) 



By substituting (fTUj) into (jS) we obtain the integral equations 



or 



a(t) = a + i^(3 f e 2i i sin( -^ds - ^ f e 2i - sm( " s) ( f e-^^VaOcfol tfa, (13) 
2 Jo 4 Jo Uo J 



a(t) = a + z— [' e 2i ^ sin ^b(s)ds, 
2 Jo 

6(t) =/5 + i— e- 2i * sin ^a(s)ds, 
2 Jo 



A 2 /•* 0;i . 



(11) 

(12) 



b(t) =(3 + i—a f e~ 2i i sin ^ds - ^ f e~ 2 ^ sm ^ 
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b(x)dx \ ds, (14) 



where a, j3 are integral constants. Therefore in the lowest order with respect to A, in 
([TO") is given by 



a + i^PJoe 2i ^ sin ^ds 



\ p + i^a^e- 2l ^ sm ^ds 
From here we consider the physical condition, ib in ([§]). therefore <j) must be normalized 
( = 1) in the lowest order, so we can set \a\ 2 + \j3\ 2 = 1. From @ the approximate 
solutions that we are looking for are 



V(0 



1 1 



V2 I ! 



V 



^i-2-sin(uit) 



+ oiJoV 



-2i-2-sin(us) 



(is | 



(15) 



where \a\ 2 + \f3\ 2 = 1. 

Finally we make a comment : If we write if) above as 

-(*)• 

then we have 

\H 2 = \ (l-cos(%m{ujt)j) (16) 

under (B) if a = (3 = e ie /^2. 

Our solutions (fT5j) are more or less well-known. However our (simple) method will 
become a starting point to find more complicated approximate solutions or to generalize 
the model from 2-levels to n-levels, [T4] . 

The author would like to Marco Frasca for several useful comments. 
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